Reversibility of Red blood Cell deformation. 
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The ability of cells to undergo reversible shape changes is often crucial to their survival. For Red 
Blood Cells (RBCs), irreversible alteration of the cell shape and flexibility often causes anemia. 
Here we show theoretically that RBCs may react irreversibly to mechanical perturbations because 
of tensile stress in their cytoskeleton. The transient polymerization of protein fibers inside the cell 
seen in sickle cell anemia or a transient external force can trigger the formation of a cytoskeleton-free 
membrane protrusion of /im dimensions. The complex relaxation kinetics of the cell shape is shown 
to be responsible for selecting the final state once the perturbation is removed, thereby controlling 
the reversibility of the deformation. In some case, tubular protrusion are expected to relax via a 
peculiar "pearling instability" . 



Red Blood Cells (RBCs) have been extensively stud- 
ied by physicists as a relatively simple example of bio- 
logical cells [H [2]. Their mechanical properties reflect 
the structure of the cell interface, including the (fluid) 
plasma membrane (PM) and the cytoskeleton (CSK), 
a two-dimensional network of flexible spectrin fllaments 
connected to the membrane through node complexes [Tj. 
This composite mechanics is responsible for a remark- 
able variety of equilibrium shape (stomatocyte-discocyte- 
echinocyte) [2J|3]. RBCs also exhibit a complex dynami- 
cal response to mechanical stress, characterized by visco- 
elasticity, and plastic deformation at high strain [5] . 

Some RBCs undergo major deformation during the 
respiratory cycle, and it appears important to under- 
stand the extent to which such large deformation might 
be reversible. In sickle-cell anemia, a mutation in the 
hemoglobin (Hg) gene leads to the formation of Hg fibers 
inside the cell and results in sickle shaped cells The 
fibers depolymerize in the lung because of oxygen intake, 
but polymerize again as oxygen is released. The cells 
loose their flexibility and may obstruct small blood cap- 
illaries, with serious medical consequences [7]. Unpub- 
lished observations suggest that in some case, rapid fiber 
depolymerization produces a large membrane bleb that 
fails to reincorporate the cell [B[ . A (very) long membrane 
tether extracted from a RBC by a external force {e.g ap- 
plied by optical tweezers) may also sometimes fail to re- 
tract into the cell when the force is switched off [S] ■ Such 
events, if they occurred in-vivo, could lead to repeated 
loss of cell membrane area and considerably shorten the 
lifetime of RBCs. 

Our goal is to determine the conditions under which 
a transient mechanical perturbation can give rise to ir- 
reversible cellular modifications. We study the model 
sketched Figjl] of a highly deformed RBC presenting a 
tubular protrusion and investigate the cell's relaxation af- 
ter the force that created the protrusion is removed. Such 
deformation can be generated by the polymerization of 
long fibers [TDl , or the application of a localized mechan- 



ical force on the cell membrane [TT]. We first show that 
the tight mechanical coupling between the CSK and the 
cell membrane can generate mechanical frustration and 
the appearance of several meta-stable cell shapes. We 
then show that the kinetics of shape relaxation strongly 
influences the relaxed shape. Finally, we briefly discuss 
more complex relaxation routes, including the peculiar 
pearling of a long and thin protrusion. 

Metastable shapes of a RBC The interplay be- 
tween CSK elasticity and the semi-permeable nature of 
the PM (permeable to water, but not to Hg or other 
larger molecules) is known to give rise to a rich phase 
diagram of equilibrium shapes, which can be explored by 
varying external parameters such as osmolarity or tem- 
perature P]. Here, we are not concerned by this kind 
of equilibrium shape transition, but wish to study the 
existence of alternative metastable shapes, where some 




FIG. 1. A strongly deformed Red Blood Cell (top) con- 
tains a long, cytoskeleton-free membrane tube filled with 
hemoglobin, created by an internal force (the polymerization 
of an hemoglobin fiber) or external force (the action of opti- 
cal tweezers). Upon force removal (/ — >■ 0) the cell may relax 
to its initial shape if the tube completely retracts (A), or to 
a different state where some membrane area remains outside 
the cell body and forms a spherical bulge (B). In the latter 
case, the tube might exhibit pearling (C) during relaxation. 
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membrane separates from the cell body, or fail to rein- 
corporate it when extracted. For this purpose, we assume 
that the reference state for the cell is a simple sphere of 
volume V, and that the total cell volume remains con- 
stant during shape transformation (we have check that 
relaxing this constraint does not alter our results). The 
elastic energy of the reference state includes contribution 
for the elastic energy of the CSK (assumed quadratic: 
~ Ka{S — 5'o)^/2S'o, with a stiffness Kc and a reference 
area 5*0, and where S = {^6ny/^V^/^ is the CSK area). 
It also includes a contribution for the plasma membrane, 
which we simply write = aS, assuming the membrane 
tension a to be constant for simplicity. We have checked 
that including a more accurate membrane elasticity jl2j 
does not alter our conclusions. 

The CSK is known to contract when the entire RBC 
membrane is removed [13] . which suggests that it is 
stretched by its attachment to the cell membrane {S > 
So). The formation of a membrane bulge detached from 
the CSK would thus reduce the CSK stretching energy, at 
the expense of the membrane deformation energy [21 [H] . 
The energy of a protrusion of area S and volume V con- 
nected to a spherical body of volume V — V is 



1^ ((36^)i/3(V-F)2/3-5o) 
2""^ So 



(1) 



The energy of a cell with a spherical protrusion Egph — 
i'ei|s=(367rV2)i/3, presented Figj2| shows the existence of 
two (meta) stable cell shapes if the CSK is sufficiently 
stiff or sufficiently stretched by the cell membrane in the 
reference state. The latter effect is characterized by the 
CSK prestress parameter: = 1— (Sq/S)'^^'^ (the excess 
cell volume compared to the optimal volume for vanish- 
ing prestress). The location of the critical line is ap- 
proximately: eo = 2^/^{a/KcY^^ — /2Kc- In practice, 
a/K^ ~ 0.1 (ct ~ 10-6j/m2 and ~ 10-5j/m2) [2]) 
and extracellular membrane protrusion can be stable at 
fairly weak prestress (eo — 25%). The shape transition is 
second order, characterized by a large (~ SOOfc^T) energy 
barrier. It does not occur spontaneously near the critical 
line, but can be triggered, e.g by an external force. 

Relaxation kinetics We assume that a localized force 
has created a long tubular protrusion of length L and ra- 
dius Ti (with L ^ ri). The force is switched off at time 
i = 0, and the perturbation relaxes. We assume for now 
that the protrusion shape relaxes smoothly toward to- 
ward a spherical bulge (without pearling), and can be 
described by two parameters (its volume V and area S) 
continuously evolving from the initial values of a thin 
tube [Vi ^ r^L and Si ^ ViL) to those of a spherical pro- 
trusion (y ^ S^^^). Cytosol volume and membrane area 
are transferred between the cell body and the protrusion 
during relaxation, leading to energy dissipation. If dissi- 
pation is dominated by cytosol hydrodynamics, volume 
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FIG. 2. a) Static phase diagram in the parameter space K^/o 
(elastic ratio) and eo (the CSK prestress), showing the region 
where a (meta) stable protrusion can exist (shaded grey). The 
dashed line is the approximation given in the text, b) Energy 
of a RBC with a spherical protrusion as a function of the 
relative protrusion volume v, for different eo {Kc/cr = 10). 
An energy barrier exist &t v = v* beyond a critical prestress. 



exchange is slow and the protrusion area decreases un- 
der almost constant volume to form a small sphere. If 
dissipation is dominated by membrane flow, the protru- 
sion volume increases under almost constant area to form 
a large sphere. When the protrusion becomes a sphere, 
the system evolves in the energy landscape shown Figj2] 
which may exhibits an energy barrier for a protrusion 
volume V* {v = V/V is the normalized protrusion vol- 
ume). The initial relaxation dynamics thus determines 
on which side of the barrier the system falls, and whether 
the final state is an intact cell (w < w*), or a cell with a 
"permanent" spherical bleb {v > v*). 

In order to study this situation more quantitatively, we 
write the balance of generalized elastic and dissipative 
forces using a Lagrangian description [T21 [H] . 



dE 



dV 



d{S,V} d{S,V) 



= 



(2) 



including elastic forces derived from the system's energy 
E and dissipative forces derived from the energy dissipa- 
tion (per unit time) V . The energy includes the CSK and 
membrane tension energy Ec\, Eq.Q, and a contribution 
from the membrane bending rigidity k, approximated by 
the bending energy of a tubular protrusion (the only limit 
where it is relevant): E^, — n/^S'^ /V'^. 

The energy dissipation functions for membrane and cy- 
tosol flows are derived in detail in the Supplementary In- 
formation (S.I.). The former is dominated by the mem- 
brane flow over the (immobile) cytoskeleton inside the 
cell body as the protrusion area varies. It is characterized 
by a friction parameter a\ ~ [(10^ — 10^^)]Pa • s/m which 
dependent logarithmically on the system's size [T7]. The 
latter is dominated by the cytosol flow through the neck 
connecting the body and the protrusion, assumed to be of 
small dimension of order the CSK mesh-size r„ ^ lOOnm 
[18]. It is characterized by a viscosity rj « lO^^Pa/m 



Vv = 



7rr; 
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Inserting E = E^x + E^, Eq.Q, and V = Vv +Vs, 
Eq.(|3|, into Eq.([2]) yields two coupled dynamical equa- 
tions for S and V. Using normalized variables and 
omitting various numerical prefactors and less important 
terms for readability {x is the normalized form of x, all 
details and the complete equations can be found in the 
S.I.), these equations read: 



ill 



-0 
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where s and v are the area and volume of the protrusion 
(normalized by those of the intact cell) and sq = (1 — 
eo)^'''^ is the optimal CSK area. The important elastic 
parameter in Eq.Q is the ratio of CSK to membrane 
elasticity Kc = Kc/cr (~ 10). The bending rigidity k = 
K/{aS) only intervenes at the very early stage of tubule 
relaxation. The important dynamical parameters are the 
typical exchange times for area (r^ ~ aiS /a) and volume 
(t„ ~ ■qS'^ / {(jr^)) . The latter depends quite strongly on 
the size r„ of the neck connecting the protrusion and the 
cell body. Interpreting Eq.Q is rather straightforward; 
once the force that caused the protrusion disappears, area 
is drawn toward the cell body (s < 0) by membrane 
tension and to some extent by membrane bending stress, 
volume is pushed out of the cell (u > 0) by the CSK (if 
it is extended - w < eo)j by the membrane tension, and 
to some extent by the bending stress in the protrusion. 

The dynamical equations Eq.Q are highly non- linear 
and must be solved numerically. A spherical protrusion 
is formed once s{t) = u(t)^/'^ (we call the corresponding 
volume Vy^ , at which point v and s are no longer indepen- 
dent, and the protrusion evolves as a sphere described by 
a single dynamical equation (given in the S.I). The shape 
of a cell with a spherical protrusion evolves in the energy 
landscape of Figj2| which possesses two minima and a 
barrier &i v — v* for large enough prestress. If > w*, 
the deformation is irreversible and a permanent protru- 
sion remains after relaxation. FigjSja illustrates the im- 
pact of the cell dynamics on the critical tether volume be- 
yond which cell shape change is irreversible. As discussed 
earlier, the protrusion is most likely to relax toward a 
large bulge with Vy^ > v* if cytoplasm exchange with the 
cell body is fast compared to membrane exchange. For 
very fast cytosol exchange (t^/ts — )• 0), the protrusion 
evolves with a constant area, and the irreversible shape 
change can be triggered by extracting a protrusion of vol- 
ume Vi > V* 1% of the cell volume. This can be done 
by extracting a lOOnm-radius membrane tether of length 
L ^ lOO^m [201 (with V = lOO^m^), a large value, but 
within experimental range. Slower cytosol to membrane 
dynamical ratio (increased t^/ts) renders cell deforma- 
tion more reversible. For RBC, on expect Ty ^ Tg and 
the relaxation time to a spherical protrusion taking up 
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FIG. 3. Critical initial protrusion volume beyond which a 
transient protrusion fails to retract inside the cell body; a) for 
a protrusion created by an external force, as a function of the 
ratio of the volume to area exchange times (with Kc/(J = 10, 
eo ~ 1/4), and b) for a protrusion created by a bundle of 
protein fibers (proteins of size 4nm), with t^/ts — 15), as a 
function of the number of fibers in the bundle, for different 
values of the protein diffusion coefficient. 



1% of the total cell area is of order 10 seconds (sec S.I). 

The case of a depolymerising fiber In the phys- 
iologically relevant case of a protrusion formed by the 
polymerization of a protein fiber (an Hg fiber in sickle cell 
anemia), which is then rapidly depolymerize at t = 0, the 
Hg concentration and the osmotic pressure are initially 
be vastly larger in the protrusion than in the cell body. 
There is a strong driving force to increase the protrusion 
volume, and irreversible shape changes are much more 
likely. Lets consider a protrusion containing N of the to- 
tal Af monomers of the cell. Assuming that the density 
difference between the two compartments vary over the 
size r„ of the neck connecting the protrusion to the cell 
body (r„ <C L), Fick's law predicts an equilibration ki- 
netics of the form: N ~ Drn ( jr — ) ■ The normalized 
kinetic equation for n e 



N/J\f reads 
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The system's free energy appearing in Eq.([2]) must also 
include the entropic contribution of both compartments: 
E = Ed + Ek. + En with (using the ideal gas law) 



/ y2 1 — Tl 

En = k^TJ^ „log - (1 - n) log 

\ V 1 — V 



est. 



(6) 



Immediately after fiber depolymerization, the osmotic 
pressure in the protrusion is very large, leading to a fast 
inflation of the tether (by volume transfer from the cell 
body with our assumption of constant total volume V). 
This additional force thus increases the likelihood of end- 
ing with a permanent protrusion, and the more so if 
the protrusion volume is large, ii) the protein density in 
the protrusion is large, and in) protein diffusion through 
the connecting neck is slow: tq/ts ^ 1. In sickled RBC, 
protrusions are formed by Hg fibers bundled together. 
Figj3]b shows the minimal volume a protrusion must have 
in order not to retract after fiber depolymerization as a 
function of the number of fibers in the bundle, for differ- 
ent values of the protein diffusion coefficient. 
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Note that water can permeate through hpid mem- 
branes [2l], and the more so through the membrane of 
cells, which contain water channels [22]. Water per- 
meation can be readily included in Vy in Eq.([3]) (with 
'Pporm ^ UpV^/S; Up a permeation coefficient), and does 
not qualitatively modify our conclusions. It is however 
likely to have a strong quantitative effect if the osmotic 
pressure difference between the protrusion and the ex- 
tracellular medium is large. Increase of the protrusion's 
volume by osmotic water influx should strongly promote 
irreversible shapes change for fiber-induced protrusions 
or under hypotonic conditions. 

Pearling of a RBC tether The previous analysis 
showed that if a sufficiently long membrane tether is ex- 
tracted from the cell by a transient force, elastic stress in 
the CSK can prevent the membrane to reintegrate the cell 
when the force is switched off, eventually leading to the 
formation of a "permanent" spherical protrusion. The 
protrusion may however not relax directly to a spherical 
shape. Experiments using optical tweezers have shown 
evidence of a pearling instability, during which the tube 
grows as a necklace of slowly inflating bubbles [HI E] . 

The most favorable case for irreversibility is when re- 
laxation occurs with little area exchange {S ~ const.). If 
the areas of both leaflets of the protrusion membrane (the 
PM is a lipid bilayer) are also constant, any variation of 
the initial curvature Ci — l/r^ of the protrusion creates 
differential stretching in the two leaflets, with an energy 
Ebl = KsS/2{hC — hCi)'^, where C is the local mem- 
brane curvature, h (~ 5nm) the membrane thickness and 
Kg (~ O.IJ/m^) the bilayer stretching modulus. The dif- 
ferential stretching effect is equivalent to the membrane 
having a spontaneous curvature Ci which hinders the re- 
laxation toward a spherical protrusion. 

Assuming that the protrusion (initially a cylinder of 
radius , length L) relaxes with time as a string of closed- 
packed bubbles of radius r{t) — 2/C{t), the driving force 
for protrusion growth is fy = —dEsL/dV + Pcff- Here, 
Pcff includes contribution from the CSK tensile stress, 
and the osmotic pressure difference with the extracellular 
medium is water permeation is present. The curvature 
and volume of a string of spherical bubbles of constant 
total area are linked by ^ '-^ L/{CiC). The driving force 
for protrusion growth thus reads: 



(7) 



If Pcff > Kgh^Cf, differential stretching is unimportant, 
and the protrusion grows steadily toward a sphere. If 
not, the force vanishes for a string of bubble of size of 
order r^. With ~ lOOnm, pearling is expected below 
a pressure of order IkPa, or an osmotic imbalance of or- 
der ImMol. Protrusions generated by protein fibers or 
in hypotonic environment should relax directly toward a 
spherical shape, but a tether drawn from a cell close to 
iso-osmolarity should relax through pearling (FigjlJC), 



provided it is long enough for the assumption of con- 
stant area to be reasonable. This phenomenon is indeed 
observed experimentally [9l [11] . In the absence of any 
membrane exchange with the cell body, relaxation dur- 
ing pearling is controlled by inter-leaflet lipid exchange 
(flip-flop) and can be slow. It should be noted that RBCs, 
like in most other cells, maintain a preferred lipid asym- 
metry lipid between the two leaflets by ATP-dependent 
flipases [2] . Whether these proteins are mechanosensitive 
and would naturally actively reduce the bilayer stretch- 
ing asymmetry is unclear. If this is so, the kinetics of 
relaxation by pearling could be strongly ATP-dependent. 

R. Briehl, P. Nassoy, and M.S. Turner are gratefully 
acknowledged for stimulating discussion and sharing un- 
published data. 
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